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HANDLING QUALITIES OF LARGE FLEaIBLE 
CONTROL-CONFIGURED AIRCRAFT 


Introduction 

This Is the first semi-annual status report on Grant No. NSG 4018. 

The work began in January, 1979 with the appointment of Mr. Supat Poopaka 
as a one-half time Graduate Research Assistant. He is a Ph.O. student 
in the School of Mechanical and Aerospace Engineering. Mr. Poopaka has 
extensive background and capability in the analytical methods of Modern 
and Optimal Control Theory, but had no course work or background in 
aircraft flight dynamics and aeroelasticity. As a result, much of this 
first reporting period has been a learning process for him. He completed 
an Aircraft Stability and Control course taught by Dr. Swaim in the 
spring semester ond in addition has rapidly become knowledgeable in 
handling qualities, pilot modeling and aeroelasticity through Intensive 
self-study. 

Discussion of Progress 

As described on pages 18 and 19 of the proposal for this work (Ref. 1), 
our approach to an analytical study of flexible airplane longitudinal 
handling qualities is to parametrically vary the natural frequencies of two 
symmetric elastic modes to Induce mode Interactions with the rigid body 
dynamics. Since the structure of the pilot model is unknown tor such 
dynamic interactions, the optimal control pilot modeling method is being 
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applled (Ref. 2) and used In conjunction with a pilot rating method 
(Ref. 3). A pole placement algorithm Is also used to ?i1nta1n rigid 
body dynamics at acceptable values on short-period and phugold fre- 
quencies and damping ratios. This should ensure that the pilot ratings 
are based on the relative amplitudes of rigid and elastic pitch 
responses and not on poor rigid body dynamics. 

Figure 1 Is a block diagram depiction of how the optimal pilot 
model Is structured and fits Into the aircraft and display dynamics 
blocks. The tracking task Is to maintain a reference rigid pitch angle 

0 ■ 0 In the presence of random turbulence as a disturbance Input. Table 

1 shows the model parameters; Table 2, some response equations; Table 3, 
the model parameter values; and Table 4, the matrix equations yielding 
the quadratic optimal control solution. 

Our Intent has been to use the optimal pilot model results to 
establish separation boundaries delineating when the pilot can discern 
rigid pitch 0 'from the total pitch angle 0^ as viewed on a flight director 
display or on the outside horizon, where 0^ Is given by 

0^(t) * 0(t) - .025Ci(t) - .029Cj(t) 

and includes the pitch contributions from two low frequency elastic modes. 
The aircraft dynamics being used are basically the B-1 airplane at a sea 
level, 949 ft/sec flight condition. 

We have completed the computer programs required for this effort and 
have preliminary results for three combinations of first and second elastic 




Figure 1. Optimal Control Model of the Pilot 
















TABLE 1. DESCRIPTION OF THE OPTIMAL CONTROL MODEL PARAMETERS 


MODEL PARAMETERS 


X Aircraft motion var1ab'’es 


u Aircraft control variables 


A Aircraft dynamics (stability deriva- 
tives with Inertial and elastic 
coupling) 


B Aircraft control effects (sensitivity 
to control and throttle deflections) 


C Motion variable display transforma- 
tion 


RELATION TO PILOT PERFORMANCE 


Pilot must observe this well enough to 
comnand aircraft and provide stability 


Pilot must use this to command aircraft and 
to provide stability 


Aircraft must be stable enough for pilot 
to control 


Aircraft must respond to the commands In 
such a way that the pilot can understand 


Displayed motion variables must be 
sufficient for command stabilization 


D Control variable display transfor- 
mation 


This Infers control observation 



Dynamic model "learned" by 
the pilot. Including neuro- 
muscular lags 


The better the pilot's knowledge of the 
aircraft and his own capabilities, the 
better he can cope with noisy measurements 


L Control gains which transform pilot's 
estimates of aircraft motions to 
control actions 


Pilot attempts to tradeoff aircraft motions 
with available control power 


V Covariance matrix of observation 
^ noise 

V Covariance matrix of neuromotor noise 
m 


Large values decrease accuracy of observa- 
tions 


Large values indicate pilot is having 
difficulties controlling aircraft 





TABLE 2. OPTIMAL CONTROL MODEL OF PILOT RESPONSE 


Aircraft Short Period Dynamics 
Augmentation System 
Mission Phase 

Vestibular Afferent Dynamics 

Controlled Element Dynamics 
X ■ Ax + Bu 

Displays 

(Visual and Vestibular) 

Measurement Vector 
y ■ Cx * Du 

Mission Requirements 
VMC/ IMC Cues 

Cost Functional Weightings 
J » E -{/{y^Qyy + u^QyU + u^QpU)d^ 

Neuromuscular/Manipulator 

Dynamics 

First Order Lag 

T„ii + u ■ m 
n m 

Disturbance Environment 

(Turbulence, Shear, Guidance 
Error) 

Shaping Filters in State Equation 
X » Ax + Bu + Ew 

Cov(w) ■ I 

Work Load Expectations 
Attention Allocation 

MIN J, Ef * f , f ^ 0 
f c.^ c c^ 


VMC = Visual Meteorological Condition 
IMC = Instrument Meteorological Condition 









TABLE 3. MODEL PARAMETERS SELECTION 
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( 5 *)* 
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_J 
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T ■ .2 s (a large displacement manipulator, high force gradient 
^ device such as elevator control on a large transport 

aircraft) 


Qp Is chosen to provide T|^ * .2 s 


Controller Time Delay t * .2 s (Typical value for human opera or 

Observation Thresholds (lOX Full Scale Value) 

- 2® , TH- - TH = 2®/s 
y y y ^ 


Observation Noise Ratios 
Additive Motor Noise 
Attention Allocation 


-20db Is a typical value 
-25db Is a typical value 



The Attention Allocation can be optimized with respect to J, 1.e., 
MIN J , but for the case being studied J is not sensitive to f ; 

therefore, the fixed value of f will be used throughout the study. 



TABLE 4. PILOT MODEL SOLUTION EQUATIONS 


Pilot Control 
Control Gain Matrix 


Pilot Kalman Esti- 
mator 


Pilot Predictor 


m(t) ■ -Lx(t) 


- (3„-‘ bIk. , A. . B. 

Co - {c|d) 








x(t - t) 
G(t - t) 




x(t - t) 
u(t - t) 


+ B,m(t - t) 


+ J:cIv:‘iyJt) - c. 


- 

|; : :]; 


A,E + EAj + W - Ecjv "‘CoE = 0 


A B 


0 -L 


. B, 


- 1 


'N 


, W 


EE' 




x(t) 

D(t) 


C(t-) + e 


M.T 


X(t - T) 
S(t - T) 


C(t - T) 


c(t) « A,C(t) + B,m(t) 


Pilot Transfer 
Function 


u(S) » -{T^S + I}‘Hf + I}"\L|0}e^^» ‘ - A)'*EC,^V^‘*y(S) 

F = {L|0)|C^^» ■ ^^^^|[SI - A)"‘ - (SI - A,)’] + (SI - AJ'I^B, 
A •= Aj - i:Co^y-‘Co 




Table 4 Continued 


Covariance of Pilot 
Predicted State 


Pilot Observation 
Noise Covariance 


Pilot Neuromotor 
Noise Covariance 



o. 




'do 


A • Ao - 

t[yj(t)] « (CJCCo^)^^ 1-1 Oy 

Efn.|(t)} r (T-'(L|0)X(l|0lV)„ . , . , 


- Vy6(t -o) 
Py 

(Vy),1 ■ 7,-’- 8| ,i 

a. » o^/ti(o^) 
o, = [E(y;(t)]'“ 


N (o^) * Describing Function Gain of Threshold 
= erfc (a^/o-VT) 

e’ - Full Attention Noise Ratio 

>1 

f^ » Attention Allocation for 


,1 = 1 ,..., 

= Motor Noise Ratio 
"'i 



TABLE 5. CASES CONSIDERED 




(i) , CO 
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in- vacuum elastic mode undamped natural frequencies 


Cl' , b] 

>e 2 e 


le 2e 

aerodynamical ly and 
ratios 

inertial ly coupled elastic mode 

undamped natural 

frequencies 

and damping 


Ul 

Sp 

Sp 

■^P 


bJ 

le 

Se 

b> 

26 

^26 
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Rigid Body 
(Bare Airframe) 

3.02 

.54 

.050 

.080 

— 

— 

— 

-- 

— 

- 

Case 1 

(Bare Airframe) 

2.83 

.53 

.072 

.021 

13.24 

.050 

21.40 

.021 

13.32 

21 

Case 2 

(Pole Placement) 

3.00 

.50 

.050 

.080 

4.00 

.047 

22.00 

.009 

3.00 

21 

Case 3 

(Pole Placement) 

3.00 

.50 

.050 

.080 

4.00 

.047 

5.00 

.009 

3.00 

4, 


